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Dear Ms. Blanchard: 

During the period of this grant two main tasks were performed: a determination of a selection 
criterion for starburst galaxies most likely to emit X-rays, and performance of a pilot study of the 
X-ray emission from nine such systems. Starburst galaxies may be expected to emit flat-spectrum 
X-ray at energies above 10 keV resulting from the various remnants of the short-lived massive stars 
which characterize the starburst. 

The investigation to determine the optimum sample resulted in a change from an X-ray selected 
(HEAO-2) sample to infrared selection based on the IRAS catalogue. A much broader sample was 
thereby available for study, and selection could be limited to only the nearest objects and still obtain 
a reasonably large sample. A sample of 99 of the brightest infrared starburst galaxies was settled on 
for the X-ray survey. 

For a set of practical size, this was then reduced to a subset of 53, based on luminosity and 
nearness. X-ray emission from these objects was individually measured from the UCSD HEAO-1 
all-sky survey in four energy bands between 13 keV to 160 keV . This data base consists of about 20 
optical disk volumes. Net significance for the result was roughly two sigma, and a very hard spectral 
shape is indicated for the net spectrum of the surveyed galaxies. With the possibility of detection 
of the class, it was then felt worthwhile to examine fluxes from these sources in other archival data. 
This was performed with the HEAO-1 A2 data and the HEAO-2 (EINSTEIN) main archive and slew 
survey. Positive results were also obtained for the sample, but again at weak significance. i t h three 

independent measures of weak X-ray fluxes from nearby starburst galaxies, we wrote a letter to the 
Astrophysical Journal (enclosed) discussing these results and their likely significance, in particular, for 
the contribution to the cosmic diffuse x-ray background, perhaps as much as 25%. A follow-up paper 
with more thorough discussion of the results and methods has been submitted to the Astrophysical 
Journal but not yet accepted for publication. We expect such acceptance shortly. 
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Principal Investigator 
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ABSTRACT 

We report the results of an investigation of X-ray emission from a sample of 53 /IMS-selected candidate 
starburst galaxies. Superposed soft and hard X-ray emission from these galaxies in the Einstein IPC and 
HE AO I A-2 and A-4 energy bands, which span 0.5 to 160 keV, is detected at the 99.6% confidence level, 
after allowing for confusion noise in the HE AO 1 data. Above 15 keV the confidence level is 97%. The com- 
bined spectrum is flat, with a (photon) power-law index of 1.0 ± 0.3. The contribution of the population of 
sources represented by this sample to the 3-50 keV residual cosmic X-ray background is estimated to be at 
least ~4% assuming no evolution. Moderate evolution, for which there is some observational evidence, 
increases this fractional contribution to ~26%. 

Subject headings: galaxies: X-rays — X-rays: sources 


1. INTRODUCTION 

Discrete sources are known to contribute significantly to the 
cosmic X-ray background (CXB), but probably not all classes 
of contributing sources have been identified. Neither do we 
know yet whether, when, and how any appreciable diffuse 
emission has been produced. Interest has traditionally been 
focused on the isotropic X-ray background spectrum in the 
energy range 3-50 keV, because of its accessibility to observa- 
tion and also because its strength and hard spectrum, unlike 
any class of known sources, pose important theoretical chal- 
lenges. The most comprehensive observational information on 
the CXB had been obtained by measurements made with the 
A-2 (Marshall et al. 1980) and A-4 (Rothschild et al. 1983; 
Gruber et al. 1984) experiments aboard the HEAO 1 satellite. 

It is thought that most of the discrete source contribution is 
due to active galactic nuclei (AGNs) whose mean power-law 
spectrum has an energy index a AGN - 0.7. Leiter & Boldt 
(1982) estimated that about 30% of the CXB flux in the band 
3-10 keV is due to AGNs. Models in which all of the CXB is 
accounted for by emission from AGNs with suitable spectral 
properties have been proposed (Schwartz & Tucker 1988; 
Fabian et al. 1990) A recent Ginga detection of a flat com- 
ponent in the spectrum of Seyfert galaxies at energies > 10 
keV (Matsuoka et al. 1990) lends some support for a larger 
AGN contribution to the CXB (Morisawa et al. 1990). Clusters 
of galaxies contribute a few percent to the flux in the above 
energy band, and probably even less at higher energies 
(Rephaeli, Gruber, & Rothschild 1987; Rephaeli & Gruber 
1988). After subtraction of the extrapolated contributions of 
AGNs and clusters, the residual CXB 3-50 keV spectrum is 
found to be surprisingly flat|(see the review by Boldt 1987). A 
basic goal is the interpretation of the flat residual CXB spec- 
trum. 

Starburst galaxies (SBGs) constitute a major class of bright 
far-infrared (FIR) sources (Soifer et al. 1986), whose enhanced 
output is probably due to warm dust heated by hot, massive 
stars. Hard X-ray binaries and supernova remnants (SNRs) 
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should also be more abundant during the starburst phase. Cor- 
respondingly, it has been suggested (Bookbinder et al. 1980; 
De Zotti 1987; Weedman 1987) that SBGs emit X-rays more 
intensely than normal galaxies and therefore may contribute 
appreciably to the CXB. 

To estimate the contribution of SBGs to the CXB, Weedman 
(1987) used his 60 ^m SBG luminosity function and a flux ratio 
of /(60 ptm)/f( 2 keV) = 2.7 x 10 8 based on Einstein measure- 
ments (Fabbiano, Feigelson, & Zamorani 1982), finding that 
SBG contribute at least 13% to the CXB flux at 2 keV. This 
estimate for the fractional contribution of SBGs to the CXB at 
2 keV is interesting and points to the possibility that a major 
new discrete component of the CXB has been identified. Obvi- 
ously, much has to be done in order to test and quantify this 
hypothesis: major unknowns are the characteristic X spectra 
and evolution of SBGs. 

In this Letter we report the results of a first study of the 
largest sample yet of SBG spectra in the combined Einstein 
IPC, and HEAO 1 A-2 and A-4 energy bands, spanning the 
range 0.5-160 keV. We use H 0 - 50 km s" 1 Mpc 1 through- 
out. 


2. THE SAMPLE 

A selection of a SBG sample would have been a simple 
matter had we had a (e.g., flux-limited) complete sample of 
SBGs. This is not yet the case. We therefore use the partial 
knowledge of SBGs to motivate our selection. An iridicative 
result from analyses of IRAS and optical observations is that 
the //MS-selected sample of SBGs is not the same as the opti- 
cally selected one (Lawrence et al. 1986; Soifer et al. 1986). The 
latter authors find that the FIR/blue luminosity ratio, L F1R /L B , 
of galaxies in the sample increases with L RR , but is not corre- 
lated with L B . Soifer et al. (1986) suggest that L F IR and L B are 
independent, and that the increase of their ratio with L ^ is 
due to enhanced FIR emission, not to increased visual extinc- 
tion. Enhanced FIR emission indicates higher abundance of 
warm dust, whereas intense blue emission would be indicative 
of a rich population of young massive stars. If X emission in 
SBGs is mainly due to high SN rate, we predict their X emis- 
sion to be more directly correlated with FIR, rather than with 
optical emission. 
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Motivated by this argument, we select SBGs from the IRAS 
Bright Galaxy Sample of Soifer et al. (1987). A luminosity- 
limited subset of this latter sample is chosen to ensure that the 
majority of the selected galaxies are SBGs: we select only gal- 
axies with EpjK > 10* 1 Lq. Presumably, this is a reasonable 
value for the cutofF luminosity, based on analyses of the bright 
IRAS sources (B. T. Soifer 1988, private communication). 
Admittedly, this selection criterion is somewhat arbitrary. 

This choice defines a subset of 1 16 galaxies of the total 324 in 
the IRAS Bright Galaxy Sample. Only eight of the 116 galaxies 
are in the HE AO I A-l catalog, of a total of 18 in the whole 
IRAS Bright Galaxy Sample. Of the 1 16 SBGs, 19 are classified 
as AGNs, whose CXB contribution has been accounted for 
separately. Therefore, the sample reduces to 97 SBGs. Only six 
of these 97 galaxies are in the sample selected by Persic et al. 
(1989) for A-2 analysis. Thus, the two samples are essentially 
distinct. 

The 97 sources in the total sample vary significantly in dis- 
tance and in In order to first determine the mean X 
spectrum of a better defined sample of nearby SBGs we extract 
a subsample of 53 sources with recession velocities lower than 
6000 km s 1 (including the two nearby SBGs NGC 253 and 
M82) for the first stage of the investigation whose results are 
reported in this Letter. 

3. X-RAY DATA ANALYSIS AND RESULTS 

The two full sky surveys performed by the A-2 and A-4 
experiments aboard the HEAO 1 satellite (Rothschild et al. 
1979; Matteson 1978) provided remarkably homogeneous 
scanning data in the 2-175 keV energy range. The HEAO 2 
Slew Survey also covers much of the sky, although at variable 
exposure. We describe the treatment of sky, or confusion, noise 
in the HEAO 1 A-4 data, because this has a strong effect on the 
significance of the results. 

We have examined HEAO 2 data from the Imaging Pro- 
portional Counter (IPC), where available, for these sources. 
The angular resolution amounts to negligible source confu- 
sion. Six sources were found listed in IPC fields, only one of 
which was identified. For another 41 positions, flux estimates 
were obtained from the Slew Survey. There is no IPC informa- 
tion on (only) six sources. 

The HEAO 1 A-2 data analysis was carried out according 
to the standard procedure described in Della Ceca et al. (1989); 
in particular, objects contaminated by (i.e., <6° away from) 
bright A-2 sources were discarded. In addition, the very selec- 
tion procedure of our SBG sample ensures that all the objects 
are at high Galactic latitude (b ^ 20°), so they are free from 
Galactic contamination. The final A-2 sample consists of 45 
objects. Their statistical error-weighted, mean individual flux is 
0.145 ± 0.041 R1 5. 

The two Low Energy Detectors (LED) of the A-4 experiment 
aboard HEAO I were sensitive in the 13-175 keV band, with 
geometric area of 103 cm 2 each, mean efficiency of 0.7, and a 
field of view 1?43 by 20°. A given source was scanned for a live 
time of the order of 1000 s at half-year intervals. For most 
sources three complete scans were obtained. A local back- 
ground count rate was also accumulated for about 13 s or 3° of 
scan azimuth preceding and following each source transit. Net 
source minus background spectra and error estimates from 
counting statistics were then obtained by co-adding the data 
from individual transits. The pulse-height data were grouped 
into four broad channels with boundaries at 13, 25, 40, 80, and 
160 keV. 
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In addition to Poisson counting noise, measurement noise 
arises from fluctuations of the diffuse background (or equiva- 
lently, source confusion), incomplete subtraction of detector 
internal background, and detector gain variations. A measure- 
ment of the sum of all these effects was obtained by binning all 
the scanning data into 1468 full-response pixels. About 750 of 
these pixels were unconfused with any of the sources in the A-4 
catalog (Levine et al. 1984) or the Galactic plane. For each 
energy band, and rms fluctuation level for the 750 pixels was 
determined. These were about 0.01 counts s' 1 , or about 1% of 
the mean background count rate, and also surprisingly close to 
the single-pixel statistical error. These values are a factor of ~ 3 
above an extrapolation of the 2.3% per 26 deg 2 quoted (Shafer 
& Fabian 1983) for the 2-10 keV diffuse background rms; thus 
noise from detector systematics appears to dominate. Since 
these systematics can potentially be modeled, the observed 
nonstatistical noise can be regarded as an upper limit. The 
total noise levels, roughly double the statistical noise, reduce 
the statistical significances by a factor of ~2. Net significances, 
statistical plus measured systematic, are reported below; these' 
should be regarded as lower limits. 

We have weighted the Einstein, A-2, and A-4 data on the 
sampled sources in two different ways. First, in order to be able 
to characterize the X-ray properties of a typical SBG in our 
sample, we weighted the individual count rates by z 2 , to correct 
for the large scatter in distance. Doing so, and assuming photo- 
electric absorption corresponding to JV H = 10 22 cm' 2 (cf. 
Schaaf et al. 1989), our best-fit power law spectrum in the ranee 
0.5-80 keV is 

( F \ -(0 99 ± 0 . 27 ) 

m = (9-48 ± 3.61) xlO-^^J 

cm -2 s“ l keV 1 , (1) 

with a reduced x 2 of 1.02. This spectrum is shown in Figure 1. 
One way to take account of the intrinsic differences in the X 
properties of the sampled sources is to make the reasonable 
(albeit unproved) assumption that the X luminosity is pro- 
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Fig. 1. — The mean X-ray spectrum of a sample of 53 /fl/4S-selected candi- 
date SBGs with z 1 weighting. Crosses represent data from the Einstein IPC, 
HEAO I A-2, and HEAO J A-4 instruments, respectively, with increasing 
energy. The fluxes in the four A-4 channels have been combined together in the 
analysis. Errors shown are 1 o . The dashed line is the best -fit power-law spec- 
trum. Net confidences are P - 0.996 for all data, P = 0.985 for HEAO 1 data 
alone, and P = 0.963 for data above 15 keV. Observed counting rates have 
been converted to photon fluxes with spectrum-dependent average efficiencies. 
Photoelectric absorption corresponding to N H = 10 22 cm -2 has been 
assumed. 
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portional to L FIR . The mean sample spectrum, obtained by 
weighting individual fluxes according to L F1R , is then 


#E) = (13.54 ±4.21) x 10 



— ( 1.05 ± 0 . 22 ) 


cm 2 s 1 keV 1 , (2) 


with a reduced x 2 of 0.62. In both equations (1) and (2), the 
quoted uncertainties include statistical and sky noise errors. 

The above deduced mean SBG spectra are significantly 
flatter than the mean AGN spectrum, whose power-law index 
is 1.7 ± 0.1 (Boldt 1987). Directly integrating the flux in equa- 
tion (1) over the band 2-30 keV, we compute a mean X lumi- 
nosity of our fiducial SBG (looted at an rms redshift of 
<z 2 > 1/2 = 0.0133) to be L 2 _ 30 k (3.3i };f) x 10 42 ergs s'|. Pre- 
vious determinations of the X luminosity of SBGs were limited 
to lower energies. For example, Fabbiano (1988) calculated the 
luminosities of the two nearby SBGs NGC 253 and M82 to be 
Lj_ l0 ^ 2.5 x 10 40 ergs s’ 1 , and L t _ l0 ^ 6.5 x 10 40 ergs s" 1 , 
respectively, in the Einstein MPC 1-10 keV band. It is inter- 
esting that the luminosity of M82 is just what we would have 
predicted after accounting for the narrower Einstein band, and 
for the fact that the of M82 is about 6 times lower than our 
sample mean. 


4. DISCUSSION 

We seem to have found that a significant fraction of the 
sources in the IRAS Bright Galaxy Sample are also intense 
X-ray sources. This result could constitute an important ingre- 
dient to theoretical studies of the nature of these sources, and 
their (likely to be related) IR and X emission mechanisms. 
Here, we are primarily interested in the implications arising 
from the identification of the sources in our sample as SBG. 

The expectation that SBGs are strong X-ray sources is based 
on their enhanced SN activity, or on their high abundance of 
X-ray binaries. Bursts of SN activity can lead to enhanced 
X-ray emission either directly, through nonthermal emission 
from SNRs, or indirectly through thermal (T ~ 10 9 K) emission 
from SN-driven galactic winds (Bookbinder et al. 1980), or 
through Compton scattering of the FIR radiation field by SN 
acceleration of relativistic electrons. Our mean SBG spectrum 
is well fitted by a power law, so the hard-X emission could 
possibly be due to Compton scattering of relativistic electrons 
off the intense FIR radiation field (see also Schaaf et al. 1989). 

Electrons of energy ~ 1 GeV are required to boost a FIR 
photon to ~10 keV. The mean FIR energy density is 
—2 x 10" 12 ergs, cm -3 , i.e. ~5 times more intense than the 
cosmic blackbody radiation. To account for the mean X lumi- 
nosity of SBGs, the energy density in electrons has to be 
~2 x 10" 10 ergs cm -3 , i.e., ~10 4 times higher than in the 
Galaxy (roughly the SBG/Galaxy X luminosity ratio). A high 
density of electrons would also yield an enhanced radio emis- 
sion, though by a factor lower than the above ratio because the 
energy loss of electrons is mainly by Compton scattering. 
Indeed, there exists a definite correlation between the radio 
( L r ) and FIR luminosities of spiral galaxies. In the FIR lumi- 
nosity range relevant to SBGs, L R oc L F i R (Wunderlich & Klein 
1988). However, the possibility that hard-X emission in SBG is 
mainly due to nonthermal emission is merely a conjecture, 
which must be carefully assessed. 

The mean SBG spectra in equations (1) and (2) are consider- 
ably flatter than the mean AGN spectrum, a desirable feature if 
SBGs are to contribute significantly to the CXB. To compute 


this contribution, we need to know the local SBG density, 
which is not known very well. Our estimate of the density is 
based on the facts that the IRAS Bright Galaxy Sample covers 
only 14,500 deg 2 of the sky, and that there are 53 sources in our 
(analyzed) sample out to ~ 1 15 Mpc. Thus, the local density of 
SBGs is n 0 = 2.4 x 10" 5 Mpc -3 . To calculate the SBG contri- 
bution to the CXB, we use equation (5.53) of Boldt (1987) for 
the spectral intensity of sources whose mean luminosity per 
(measured) energy interval is dL/dE 0 , i.e.. 


dl 

dE 0 



Q) * 


( 3 ) 


where / is a dimensionless function which depends on a, on the 
maximum redshift z m out to which these sources are found, on 
the cosmological density parameter, Q 0 , and on a parameter 
characterizing the evolution of the sources, denoted by Q in 
Boldt (1987). Integration of our deduced mean spectrum (cf. eq. 
[1]) over the energy band [£ t , E 2 ] leads to 

I = <z 2 >n 0 /l/(£i 01 -£} 01 )keVcm“ I s“ 1 sr' 1 , (4) 

with A = 9.6 x 10" 5 and E x and E 2 are in keV. Assuming 
Q = 0 (no evolution), Qq = 1 and z = 4, which give/ =0.61, 
and using an rms redshift of <z 2 > 1/2 =0.0133, we calculated 
the integrated SBG intensity over the 3-10 and 3-50 keV 
energy bands to be / 3 _ 10 ^ 0.37 and / 3 _ 50 - 2.5 keV cm" 2 s“ 1 
sr" 1 , respectively. Density and/or luminosity evolution would 
yield higher predicted SBG intensities: e.g., an evolution at the 
level suggested by Danese et al. (1987) for radio-selected SBG 
candidates, corresponding to Q ~ 3, would raise (see eq. [5.57] 
of Boldt 1987) the above intensities to / 3 _ 10 ^ 2.4 and / 3 _ 50 a* 
16.6 keV cm" 2 s" 1 sr" 1 . 

From Figure 3.6 of Boldt (1987), we then estimate the inte- 
grated intensity of the residual CXB to be F 3 _ l0 a 20 and 
^ 3 - 5 o — 65 keV cm" 2 s" 1 sr" 1 . Correspondingly, with no evo- 
lution the SBG contribution to the residual CXB is ~2% in 
the 3-10 keV band, and ~4% in the 3-50 keV band. Moderate 
evolution at the above level, Q ~ 3, would increase these values 
to ~ 12% and ~26%, respectively. (Note that since the mean 
SBG spectrum is flatter than the residual CXB spectrum, the 
SBG contribution to the residual CXB is larger at higher ener- 
gies.) As for the contribution to the total CXB, equation (3.1) of 
Boldt (1987) yields / f 3 _ 10 s* 27 and r 3 . 50 a 88 keV cm" 2 s" 1 
sr" so the contribution of SBG to the total CXB is ~ 1% and 
— 3% with no evolution, and ~9% and ~19%, if Q — 3. It 
should be emphasized that since we have sampled only the 
brightest SBGs, the total SBG contribution to the CXB may be 
higher than estimated here. 

The above estimates confirm the earlier conjecture of 
Weedman (1987) and, more recently, of Griffiths & Padovani 
(1990) that SBGs contribute appreciably to the residual CXB. 
In particular, our estimate for the fractional contribution of 
SBGs to the CXB is in general agreement with that of Griffiths 
& Padovani (1990), whose (more phenomenological) analysis 
was, however, limited to Einstein data. Because our statistical 
study of hard-X emission from SBGs is the first of its kind, and 
given the fact that individual SBGs have not yet been detected 
at hard-X energies, our mean spectrum at energies ^20 keV 
may be regarded as tentative, pending more sensitive observa- 
tions of SBG spectra. We emphasize that since our deduced 
mean SBG luminosity is quite high, Ginga might well detect 
individual bright SBGs. 
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Assuming all the presumed 150 sources (over the full 4 n sr) are rarer; correspondingly, therefore, at higher energies the 

to be SBGs yields an upper limit to the contribution of these CXB may be expected to be less isotropic than deduced from 

sources to the CXB. Among other uncertainties, it is expected the IPC measurements. Indeed, although we have considered 

that an unknown fraction of these sources are AGNs. only the very bright SBGs, it is clear that lower luminosity 

However, that this fraction is small can be seen from an inde- sources also contribute, especially so at lower energies, perhaps 

pendent, previous estimate of the mean density of SBGs. The even relatively more than the bright sources. Since such 

density of SBGs with L FIR > 10 n L© can be deduced from e.g., sources are more numerous (as can be seen in Fig. 1 of Soifer 

Figure 1 of Soifer et al. (1986): scaling to our adopted value of et al. 1986), we do expect lower CXB granularity that can be 

H 0y we infer a value of ^2 x 10 5 Mpc 3 , in good agreement inferred from our calculated SBG surface density, 

with our estimate. An expanded discussion of all aspects of this work will be 

The angular surface density of SBGs, a a , can be calculated given in a forthcoming paper (Gruber et al. 1991). 
from the expression for the differential surface density of 
sources per unit redshift interval (which is directly related to, 

e.g., eq. [5.47] of Boldt 1987). Integrating to z m = 4, we obtain This research has been supported by NASA through grant 

o s ~ 700 SBGs deg -2 . This number is smaller by a factor -7 NAG 5-1385. Y. R. and M. P. are grateful to the University of 

than the lower limit obtained by Hamilton & Helfand (1987) California, San Diego for hospitality. Martin Elvis assisted us 

from their analysis of the isotropy of the CXB from Einstein with the Slew Survey. It is a pleasure to thank Elihu Boldt for 

IPC data. Since the IPC energy band (0.5-3.5 keV) is much many stimulating conversations. We also thank an anonymous 

below the range explored here, it is not at all clear that this referee for stimulating comments which helped us improve the 

limit should apply. It is reasonable that higher energy sources presentation of this work. 
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